Quasiparticles beyond the Fermi liquid and heavy fermion criticahty 
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We give a self-consistent theory of the scale dependent effective mass enhancement m* /m of 
quasiparticles by 3D antiferromagnetic (AFM) spin fluctuations in the presence of disorder at an 
AFM quantum critical point. The coupling of fermionic and bosonic degrees of freedom in the 
critical regime is described in terms of a critical quasiparticle theory. Using the fact that even in the 
"non- Fermi liquid" regime the quasiparticle width does not exceed the quasiparticle energy, we adopt 
relations from Fermi liquid theory to determine the dependence of the spin fluctuation spectrum 
on m*/m, from which the self energy and hence m*/m may be calculated. The self-consistent 
equation for m* /m has a strong coupling solution provided the initial value is sufficiently large. 
We argue that in YbRh2Si2 (YRS), quasi-2D AFM and/or 3D ferromagnetic Gaussian fluctuations 
existing over a wide range drive the system into the 3D strongly coupled fluctuation regime. We 
find critical exponents of the temperature dependence of the specific heat coefficient 7 oc and 
of the resistivity p{T) — p(0) -I- AT^^'*' in good agreement with experiments on YbRh2Si2 in the 
temperature range T < 0.3K. 
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Introduction. Quantum phase transitions differ from 
thermodynamic phase transitions in that temporal fluc- 
tuations are of equal importance or may be even dom- 
inant compared to spatial fluctuations. Early theories 
of quantum critical behavior, formulated in the frame- 
work of a Ginzburg-Landau- Wilson action of the order 
parameter field 0, Q , found that the effective dimen- 
sion of the corresponding 0'*-field theory is increased to 
def f — d + z where d, z are the spatial dimension of the 
fluctuations and the dynamical critical exponent, respec- 
tively. In many cases of interest de// is above the upper 
critical dimension, so that the fluctuations are effectively 
non-interacting and the theory is of the Gaussian type. 

However, several metallic compounds exhibit a quan- 
tum critical point with non-Gaussian character [3, [J (for 
a review see [^). Many of the systems for which the 
Gaussian theory is found to fail show strong electronic 
correlations, as exemplified by a strongly enhanced ef- 
fective mass ratio m* /m (heavy fermion systems). One 
might expect that in such systems it is not possible to 
integrate out the fermionic degrees of freedom in favor 
of a purely bosonic theory. Rather, in many cases the 
fermionic properties, in particular the effective quasipar- 
ticle mass, develop critical behavior, which will feed back 
into the bosonic spectrum. One is then faced with formu- 
lating a critical theory of coupled fermionic and bosonic 
degrees of freedom. Moreover, all signs indicate that such 
a theory would be in the strong-coupling regime. 

Over the past ten years several proposals dealing with 
this situation have been formulated [6|-[8| . Some of these 
scenarios have been developed enough to allow compar- 
ison with experimentally observed critical exponents, in 
particular for CcCus.gAuo.i and for YbRh2Si2 (9l.[lo|. 

In this letter, we analyze the behavior near an AFM 
QCP within a heavy-quasiparticle picture based on the 



Anderson lattice model and we refer our results to the 
heavy fermion compound YbRh2Si2 (YRS). When the 
energy scale Tk of the lattice Kondo effect and that of 
the onset of the critical regime Tcr are as wide apart as 
they are for YRS, namely Tk « 25A' and T„ « O.SAT, 
we may assume that well below T/f, the Kondo screen- 
ing and hence the heavy-quasiparticle picture are robust. 
However, the effect of AFM quantum critical fluctua- 
tions leads to a further renormalization of the quasipar- 
ticle spectrum, which generates a critical (scaling) be- 



havior of the effective mass [ll| . In contrast to other ap- 
proaches {e.g. [lo|). we address the non-Gaussian critical 
region and consider the interaction of the heavy quasi- 
particles with AFM fluctuations. We shall show that a 
self-consistent theory of a divergent quasiparticle effec- 
tive mass generated by interaction with 3D AFM fluctu- 
ations is capable of providing a sufficiently accurate ac- 
count of the critical region as seen experimentally. The 
interplay of spin fluctuations and fermionic excitations 
has also been considered in 1/N expansion by Abanov 
and Chubukov [l3 | and a renormalization group formu- 
lation has been given by Metitski and Sachdcv [ISj] . 

The above scenario depends sensitively on the detailed 
nature of spin fluctuations in a given system. For ex- 
ample, 3D AFM fluctuations do not lead to true criti- 
cal behavior, i.e. a Gaussian fluctuation theory is ap- 
plicable [l|, 13], provided the effective mass ratio is not 
too large. We argue below that in YRS one has a wide 
region of quasi-2D antiferromagnetic (or else 3D ferro- 
magnetic) fluctuations [Ol, which give rise to a substan- 
tial enhancement of the effective mass and drive the sys- 
tem into a strong coupling regime of 3D antiferromag- 
netic fluctuations. In YRS, the quantum critical point 
is accessed by tuning the magnetic field. At the criti- 
cal value of the field, the crossover from quasi-2D to 3D 
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antiferromagnetic fluctuations as a function of decreas- 
ing temperature takes place at T = Tcr- This crossover 
is clearly seen in experiment. Whereas in the quasi-2D 
AFM Gaussian fluctuation regime, marginal Fermi liquid 
(MFL) behavior 15| is observed experimentally (by this, 
we mean specific heat coefficient 7(T) cx h\{TolT), resis- 
tivity p{T) — p(0) oc T), in the critical regime the data 
are well-described by different power law dependences: 
7(r) oc T"-i, p{T) - p[Q) cx T". These results follow 
only when the strongly anisotropic (k-dependent) contri- 
bution of AFM fluctuations to the quasiparticle energy 
is smeared over the Fermi surface by sufficiently strong 
impurity scattering (elastic mean free path shorter than 
inelastic mean free path). 

Critical quasiparticle picture. Our starting point is an 
Anderson lattice model of correlated /-electrons (spin tj, 
site i; operators fi„) hybridizing with conduction elec- 
trons (ckcr), as described by the Hamiltonian 

k.O" i..(7 i 

+ Vj:{e^^'c+J^<r + h.c.). 

The single particle Green's functions G'ab(k, w), a,b — 
{c, /}, may be decomposed into a quasiparticle term 
and an incoherent contribution, G'(k, w) = zG^^ + G""^, 
where the quasiparticle weight factor z is defined by 
z^i = 1 - aReS](w)/9w|„=£;,. Here S(a;) is the /- 
electron self energy generated from the Coulomb repul- 
sion U and which we take to be k- independent. In the 
limit of large U and close to half- filling {rif <\ electrons 
per site), one expects z ^ 1, leading to a large effective 
quasiparticle mass ratio z^^ = m*/m. We assume that 
the Fermi level (at k — kp) intersects the lower band, 
(sec Fig. [ij . Then the quasiparticle energy is given by 
Ej^ ~ {m/m*)vp{k 



kp), where vp = 2{V/eki. 



the Fermi velocity of the uncorrelated hybridized band 

zImE(£'j:). That is 



and the quasiparticle width is F 
to say G"?P(k,a;) ^[uj-E- - iV]-^ . 

The condition for the quasiparticle picture to be valid 
is r < \E'^\. In the Fermi hquid regime, F = c{E'^Y ^ 
\E^ I in the limit E^ — s> 0. Here, we argue that the quasi- 
particle stability condition may be even satisfied in non- 
Fermi liquid situations. We extend the usual quasiparti- 
cle picture by allowing the parameter z = m/m* to de- 
pend on the energy scale, z = z{(jj) = 1/[1 — dT,{ijo) / duj]. 
It is important to observe that the (retarded) self energy 
is an analytic function in the upper half plane, so that the 
real and imaginary parts of any nonanalytic term (in the 
lower half plane) are locally connected. A case of partic- 
ular interest is the MFL form of the retarded self energy 
at T = E(w) = ci[cjln |a;|/wc- j(7r/2)|w|], yielding 
the ratio Tl\E'^\ = 7r/(2 In |) < 1 . Thus quasi- 

particles are still well defined in a MFL. Even in a true 
non-Fermi liquid phase with E(a;) oc — i(i|aj|)", a < 1 
so that ImS(a;) oc ReS(a;) oc and z oc (|£'jr|)^~", 




FIG. 1. Hybridized bands of the Anderson lattice model. 
Near half filling, the Fermi level Ep is near the top of the 
lower band. 



one finds T/\E^. \ = cot(|a) < 1 for i < a < 1. These 
examples make clear that even if z = at the Fermi sur- 
face, the spectral function for non-zero excitation energy 
may be peaked sharply enough to separate a quasiparticle 
contribution from the incoherent part. 

Spectrum of antiferromagnetic spin fluctuations. We 
shall assume that in the heavy-quasiparticle paramag- 
netic phase, the self energy for the single particle Green's 
function is determined by the interaction with magnetic 
fluctuations. Here, we model the imaginary part of the 
renormalized retarded dynamical spin susceptibility for 
wave vectors q near the AFM ordering wave vector Q by 



Imx(q, v) = 



(No/z){iy/v*pQ) 



[l + F{Q) + z{q-Qm? + iWv*pQr' 

(1) 

where A^o is the bare density of states at the Fermi sur- 
face, Vp = {m/m*)vp is the renormalized quasiparti- 
cle Fermi velocity, — kp^ is the microscopic AFM 
correlation length and F{Q) is a dimensionless general- 
ized Landau parameter, with F{Q) — > — 1 at the crit- 
ical point. The factor z multiplying (q — Q)^ in Eq. 
(1) arises as follows: Consider the small (q — Q) expan- 
sion of the unscreened static quasiparticle susceptibility 
X ^ (A^o/^)[l - a(q - Q)^?^]- The (q - Q)^ correction 
is governed by high-energy contributions and should not 
be renormalized, which leads to a = z = m/m*. Note 
that the physical correlation length ^ = £,o[{m* /m){l + 
F{Q))]~'^/'^ diverges when the critical point is approached 
in any direction. As described above, we shall assume 
m/m* = z ~ z(£), where £ is the relevant energy scale, 
e.g. temperature T or magnetic field H (more precisely 
the Zeeman splitting), whichever is larger. 

Quasiparticle self-energy. Now we set up a self- 
consistent determination of the quasiparticle self energy 
via the leading term in a skeleton graph expansion. Thus, 
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the quasiparticle width is given by 
rM = A2(Q) / — ^ImG"'P(k-q,c.-t.)Imx(q,i^) 
x[b{,,)-f(oj-u)], (2) 

where X{Q) = zF{Q)/Nq is the dimensionful (Landau) 
interaction vertex and f{uj),b{uj) are Fermi and Bose 
functions, which at low T confine the i/- integration to 
the interval [0, w]. This is an equation of self consistency 
since F = zlmT, appears non-linearly in the integrand. 

We consider 3D spin fluctuations and account for im- 
purity scattering by averaging the right hand side of 
Eq. (2) over the Fermi surface. Near the critical line, 
F{Q) « —1, and for Q of order kp, we find 



-1/2, 



-2^3/2 



(3) 



The dependence r(a;) oc cj-^/-^ in the case of 3D AFM fluc- 
tuations is well known. The structure of the full self en- 
ergy = r(a;)/z may be determined from this result 
as oc (iw)'^/^/z'^(a;). Here we have generalized to 

the frequency dependent z-factor. We may take a power 
law form z{uj) = buj" This enables the self-consistent de- 
termination of z{uj) = (1 — dKcT^/duj)^^: 

z-\io) = 1 + 4(fcF^o&)"'(t'FO)"'/'3(l/2 - a)wi/2-3a 

= 1 + C3ikp^orHvFQ)-'^'z-H^)^'^' (4) 

where C3 is 0(1). We explore the consequences of the 
scale dependent z. For frequencies less than the tem- 
perature, we replace a; by T . Since C3{kp£^Q)^^ « 1, 
we can say that as long as z^^{T){T /vfQY^"^ ^ 1 for 
any T , the system will be in the Gaussian fluctuation 
regime all the way down to the critical point. If how- 
ever, the initial value of z~^{T), when one enters the 3D 
AFM fluctuation regime, is sufficiently large, such that 
z~^{T){T /vpQY/'^ ^ 1, a new regime is accessed, which 
is of a strong-coupling nature. We find the characteris- 
tics of this new regime within the present approximation 
by solving the self-consistent Eq. (4), to get 



z{T) = [c^{kFS,^)-^Y/\TlvFQY/\ 



(5) 



In the case of only 3D AFM fiuctuations it is difficult 
to satisfy the strong-coupling condition unless z~^(T) is 
sufficiently large. Therefore, if on the initial approach 
to the critical point, fluctuations dominate that lead to 
a growing z~^(T) with decreasing T, the condition may 
be met at some point. The precise crossover point is de- 
termined by the crossover of these precursor fluctuations 
to the critical 3D AFM fluctuations and by the condi- 
tion above that leads to Eq. (5). As mentioned in the 
introduction, there are clear indications in the data on 
YbRh2Si2 of both quasi-2D AFM and 3D FM fluctua- 
tions. In both cases one finds z~^(T) oc ln(To/T), so 



that z~'^ grows as T — >• and is about 40 in the heavy 
Fermi liquid region of the phase diagram. 

Specific heat and electrical resistivity. Within the ap- 
proximation of neglecting the momentum dependence of 
the self-energy, the entropy density is given by S/V = 
(27Vo/T) / dujuj[-df /duj)[uj - ReS](w)]. Substituting 
the power law dependence found above, Rel](ti;) oc 
|tjj|3/4 sign(a;), we find a specific heat coefficient diverging 
in the limit T ^ 0, 



7(T) = c^N^iT/vpQr^'^ 



(6) 



A comparison of the theoretical temperature dependence 
of Eq. (6) with experiment 17| is shown in Fig. [5] To 
achieve this excellent fit, a T- independent constant spe- 
cific heat has to be added. Such a term could arise 
from very low frequency (a; <^ T) oscillators. Its mag- 
nitude represents about 0.4% of the total number of for- 
mula units. A possible source of such low frequency 
oscillators is spatially and temporally fluctuating AFM 
domains that oscillate about the preferred c-dircction. 
An anisotropic exchange interaction has been proposed 
to explain the observed electron spin resonance (7-shift 
[ist . We estimate the corresponding oscillator quantum 
as hjj K. 0.5(7 K, where cr <C 1 is the staggered magneti- 
zation per formula unit of a typical domain. 




r(K) 



FIG. 2. Specific heat: Comparison of theory, Eq. (6) and data 
of Ref. [13 at the critical magnetic field and below the critical 
temperature for quantum critical scaling. 

As mentioned above, beyond the critical regime proper, 
at temperatures T > Tcr ~ Q.iK, the data on YRS in- 
dicate the existence of Gaussian fluctuations of quasi- 
2D AFM or 3D FM character. This leads to 7(T) = 
CGNn{dldT)[T\n{Th/{T + Tfl)), where Tfl is the 
crossover temperature into the Fermi liquid regime and 
Th « 20K is a high temperature cutoff scale of the order 
of the lattice Kondo temperature. 

The electrical resistivity in the presence of impurity 
scattering may be obtained from the quasiparticle relax- 
ation rate F as /9(T) = p{0)+Cp{m/e'^n){m* /m)T . Using 
the above results in the scaling regime we find 

p{T) - p(0) = c,(m/e2„)(i-^.Q)i/4T3/4 (7) 
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Upon entering the Gaussian fluctuation regime at T > 
Tcr this fractional power law behavior crosses over into 
a linear T-dependence. A comparison of the theoretical 
temperature dependence of Eq. (7) with experiment [l^ 
at the critical magnetic field is shown in Fig. [31 
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FIG. 3. Resistivity: Comparison of theory, Eq. (7) and data 
of Ref. [l^ at the critical magnetic field and below the critical 
temperature for quantum critical scaling. 



uid theory are observed. We argue that the MFL behav- 
ior is caused by Gaussian fluctuations of quasi- 2D anti- 
ferromagnetic and/or 3D ferromagnetic character. These 
fluctuations lead to a quasiparticlc effective mass m* in- 
creasing logarithmically with temperature. By contrast, 
3D antiferromagnetic fluctuations do not lead to an in- 
creasing m*, unless m* exceeds a certain threshold. Then 
a strong coupling regime is reached. We calculated the 
critical exponents in that regime in a skeleton graph ap- 
proach. Results for specific heat, resistivity and NMR 
relaxation rate are in Eqs. (6-8). In Figs. 2,3, we com- 
pare our results with the low-T data on specific heat and 
resistivity and find reasonable agreement. 
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search unit 960 "Quantum phase transitions." Part of 
this work was carried out at Los Alamos National Labo- 
ratory (PW) and the Aspen Center for Physics (EA,PW). 



Local susceptibility. The local (g-intcgratcd) suscepti- 
bility xioc{^) determines the nuclear spin relaxation time 
Ti through the relation 



^oc[iimxiocHUo = E;;7? 



No{m*/mY 



Na{m*/m)\i/^o)/{ST:vFQeo)- 



where ^ is the physical spin correlation length already 
defined below Eq. (1) as ^ = £,o[{m* /m){l + F{Q))]-^/^ 
and we used Eq. (1). Close to, but not quite at the critical 
point, such that limT^o(l + F{Q)) > 0, we then find 

111 m 

Exactly at the critical point the temperature dependence 
of the control parameter enters, [1 -I- F{Q)] oc , and 
l/TiT oc T-'^/s-P/'2. the value of /? is not known at 
present. Eq. (8) is a prediction of our theory; unfor- 
tunately the lowest field used in the available data [2o| . 
0.15 T, is about twice the critical field of 0.06 T and is 
thus outside the critical region 

Conclusion. We have presented a critical quasiparticlc 
theory of heavy fermion compounds near a critical point. 
We focused on YbRh2Si2, but the general framework of 
our approach should be applicable to other systems. The 
starting point is the observation that quasiparticles are 
well-defined in certain classes of non-Fermi liquid states. 
The best known example is the marginal Fermi liquid, 
properties of which are observed in a certain region of 
the T — H phase diagram of YRS. Close to the critical 
point, however, much stronger deviations from Fermi liq- 
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